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Summary 

"  In  this  paper  the  fully  nonlinear  equilibrium  theory  of  homogeneous 
and  isotropic,  incompressible  elastic  solids  is  used  to  study  the  elasto- 
static  field  in  plane  strain  on  a  half-plane  deformed  by  a  concentrated 
surface  load.  Under  suitable  restrictions  on  the  form  of  the  elastic 
potential  at  severe  deformations,  it  is  shown  that,  for  materials  which 
ultimately  ’’harden*  in  simple  shear,  the  displacement  is  bounded  near  the 
point  of  application  of  the  load.  This  is  not  the  case  for  materials  which 
ultimately  "soften*  In  shear.  Estimates  of  the  true  stress  tensor  near  the 
singular  point  are  given. 
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I NTRODUCTI ON 

Under  certain  circumstances,  the  description  of  the  elastostatic 
field  furnished  by  the  classical  linearized  theory  of  elasticity  may  be 
inadequate,  even  when  the  applied  loads  are  small.  Such  breakdowns 
in  the  linearized  theory  are  ordinarily  local  in  nature  and  are  brought 
about,  for  example,  by  stress  concentrations  such  as  those  induced  by 
holes  or  cracks  in  the  interior  of  the  loaded  solid.  The  most  extreme 
examples  of  problems  of  this  kind  involve  a  singular  point  in  the 
elastostatic  field— the  tip  of  a  crack,  for  example— near  which  the  dis¬ 
placement  gradient  is  unbounded.  Since  the  basic  approximative  assump¬ 
tion  underlying  the  linear  theory  requires  that  this  gradient  be 
negligibly  small  in  comparison  with  unity,  it  is  hardly  surprising  that 
results  based  on  this  theory  may  be  in  error  near  such  a  singular  point. 

Problems  involving  large  displacement  gradients  properly  fall 

within  the  scope  of  the  finite  theory  of  elasticity.  In  recent  years 
there  have  been  several  investigations  within  the  framework  of  the 
finite  theory  of  the  local  structure  of  the  elastostatic  field  near  a 
geometrically- induced  singular  point.  Much  of  this  work  is  summarized 
in  the  review  articles  [1,2],  where  references  are  given.  In  general, 
the  analyses  of  singular  problems  reviewed  in  [1,2]  are  necessarily 
local  in  character;  they  reveal  that  the  results  from  linear  theory 
near  the  singular  point  are  Invariably  incorrect  quantitatively,  and 
In  some  Instances  may  be  qualitatively  misleading  as  well.  Since  it  is 
often  the  field  near  the  singular  point  which  is  of  primary  physical 
Interest,  analyses  based  on  finite  elasticity  are  of  considerable 
significance. 


i. 
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In  the  present  paper,  a  further  singular  problem  in  elasto- 
statics  is  considered  within  the  scope  of  the  theory  of  finite 
elasticity.  This  is  the  plane  strain  problem  of  a  con¬ 
centrated  uniform  normal  line  force  applied  to  an  elastic  body  which, 
in  the  undeformed  state,  occupies  a  half-space.  Here  the  singularity 
arises  because  of  the  character  of  the  applied  load,  rather  than  from 
the  geometry  of  the  undefomed  body,  as  is  the  case  in  most  singular 
problems  previously  treated  within  the  finite  theory  [1,2].  The 
present  analysis  aims  at  the  asymptotic  determination  of  the  displace¬ 
ments  and  stresses  near  the  point  of  application  of  the  load.  We  deal 
with  the  fully  nonlinear  equilibrium  theory  for  homogeneous,  isotropic 
incompressible  materials  that  possess  an  elastic  potential.  The 
only  restriction  on  this  potential  is  one  which;  pertains  to  its 
asymptotic*behavior  at  large  deformations;  it  is  this  regime  of 
deformation  which  dominates  the  local  field  near  the  singular  point. 
Again,  it  is  found  that  the  structure  of  the  stress  and  displacement 
fields  near  the  singular  point  differs  from  that  predicted  by  the  linear 
theory. 

The  only  previous  works  devoted  to  the  effect  of  nonlinearity  on 
the  elastostatic  field  near  the  point  of  application  of  a  concentrated 
force  are  those  of  Arutiunian  [3]  and  Atkinson  [4],  Both  of  these 
authors  retain  the  assumption  of  infinitesimal  displacement  gradients 
appropriate  to  the  linearized  theory,  but  replace  the  constitutive  law  of 
the  latter  theory  by  a  nonlinear  one, 


Section  1  contains  a  review  of  some  prerequisites  from  the  theory 
of  finite  plane  elastostatlcs  for  homogeneous,  Isotropic,  incompressible 


elastic  solids.  We  also  Introduce  In  Section  1  the  special  class  of 
such  solids  underlying  the  subsequent  analysis.  Section  2  is  devoted 
to  the  formulation,  analysis,  and  discussion  of  the  concentrated  force 
problem.  Only  the  case  of  a  tensile  force  is  treated  in  detail. 
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1.  PRELIMINARIES  FROM  PLANE  FINITE  ELASTOSTATICS 

In  this  work  we  shall  be  concerned  with  the  analysis,  within  the 
finite  theory,  of  plane  elastostatlc  fields  in  Incompressible,  homo¬ 
geneous,  and  isotropic  elastic  materials  In  the  absence  of  body  forces.^ 
Consider  an  elastic  body  which— in  the  undeformed  state— is  an 
Infinite  cylinder,  and  let  n  denote  a  plane  open  cross-section  of  this 

cylinder  perpendicular  to  its  generators.  Let  (x^.Xg)  be  the  coordi¬ 
nates  of  a  generic  point  in  H  relative  to  a  fixed  two-dimensional 
rectangular  cartesian  coordinate  system  in  the  plane  of  II  . 

A  plane  deformation  of  the  body  is  given  by  the  transformation 

ya  *  Vxrx2)  *  xa  +  ua*XrX2*  on  n  »  a=1*2»  (1.1) 

where  are  the  components  of  the  position  vector  y  of  the  particle  in 
the  deformed  body  whose  position  vector  in  the  undeformed  configuration 
is  x;  ua  are  the  components  of  the  displacement  vector  u,  all  with 
respect  to  the  rectangular  coordinate  system.  The  function  y  is 
required  to  be  twice  continuously  differentiable  on  n,  and  it  is  fur¬ 
ther  required  that  the  mapping  x  y  be  one-to-one  and  that  its 
inverse  x  have  the  same  smoothness. 

The  deformation  gradient  tensor  F  associated  with  y  has  compon¬ 
ents 

3ya 

Fa$a  9x^  ’  0-2) 

Vor  a  discussion  of  the  foundations  of  finite  elasticity  see  Gurtin  [5]. 
For  further  reading  on  plane  finite  elastostatics  of  Incompressible 
materials,  see  [6]. 
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Slnce  the  material  is  presumed  to  be  incompressible,  the  deformation 
(1.1)  must  be  locally  volume-preserving,  whence  the  Jacobian  determin¬ 
ant  of  the  mapping  must  satisfy 

J  =  det[Fae]=l  on  n  .  (1.3) 

Define  the  right  and  left  two-dimensional  Cauchy-Green  tensors  C  and  G, 
respectively,  by 

C  =  FTF  ,  3  =  FFT  .  (1.4) 

These  deformation  tensors  have  common  fundamental  scalar  invariants 
given  by 


!rtr5s  FaeFae  s  1  ’  S3*’ 

I2  *  det  C  =  J2  =  1  . 


(1.5)T 


The  invariant  I  is  found  to  obey 

I  >  2  on  II  •  (1.6) 

Moreover,  I  ■  2  if  and  only  if  F  ■  1,  where  1  is  the  two-dimensional 
unit  tensor. 

Let  t  be  the  two-dimensional  true  (Cauchy)  stress  tensor  regarded 
as  a  function  of  position  on  the  deformation  image  n*of  n.  Its  compon¬ 
ents  tag  represent  forces  per  unit  deformed  area.  If  a  is  the  asso¬ 
ciated  nominal  (Plola)  stress  tensor  field  on  n,  whose  components  aag 

^Repeated  subscripts  are  summed  over  the  range  (1,2). 
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represent  forces  per  unit  undeformed  area,  one  has 

a  =  t(FT)_1  .  (1.7) 

For  an  equilibrium  deformation  in  the  absence  of  body  forces,  it  is 
necessary  that  t  satisfy 

div  r  =  0  ,  t  =  on  H*  .  (1.8) 

It  follows  from  (1.7), (1.8)  that 


div  a  =  0 


on  IT  . 


(1.9) 


Suppose  that  r  is  a  regular  arc  in  n  which  is  mapped  onto  r*  in 
n  by  the  deformation  (1.1),  and  denote  by  n  and  n*  unit  normal  vectors 
of  r  and  r*,  respectively.  The  true  traction  vector  t  and  the  asso¬ 
ciated  nominal  traction  vector  s  are  given  by 


s  =  a  n  on  r  , 
t  !  in*  on  r*  , 

It  can  be  shown  that 


(1.10) 


s  *  0  on  r  if  and  only  if  t  *  0  on  r*  .  (1.11) 

Moreover,  (1.11)  continues  to  hold  true  for  an  arc  r  on  the  boundary  of 
n  If  the  deformation  and  nominal  stress  field  are  suitably  regular  on 
the  closure  f  of  H.  This  important  fact  allows  the  boundary  condition 
for  a  traction-free  surface  r*  in  the  deformed  body  to  be  specified  on 
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the  known  pre-image  r  of  r*  in  the  undeformed  body. 

The  mechanical  response  of  the  homogeneous,  isotropic,  incompres¬ 
sible  material  under  consideration  is  governed  by  the  strain  energy 
density  W  per  unit  undeformed  volume.  For  a  plane  deformation  of  the 
type  described  above,  W  depends  only  on  the  deformation  invariant  I: 

W  =  W(I)  •  (1.12) 

The  stress-deformation  relation  is 

'at  '  FcpFSp  -  »Sc,0  <K,3> 

where  6  „  is  the  Kronecker  delta  and  the  scalar  field  p  is  an  arbitrary 

op 

hydrostatic  pressure  whose  presence  is  necessary  because  of  the  con¬ 
straint  of  incompressibility.  Because  of  the  presence  of  p,  the 
true  stress  tensor  is  not  completely  determined  by  the  deformation  for 
an  incompressible  material .  From(1.13),  (1.7)it  follows  that 


a  a  =  2W'(I)F  .  -  peQ  e  F  on  II,  (1.14) 
aB  '  aB  By  ap  py  ’ 


provided  eag  are  the  components  of  the  two-dimensional  alternator.  In 
the  foregoing,  W  denotes  the  derivative  of  W  with  respect  to  I;  we 
assume  that  W  is  twice  continuously  differentiable  for  I  >.  2.  It  is 
further  assumed  that  W  vanishes  in  the  undeformed  state,  so  that 


W(2)  =  0  , 

and  that 


(1.15) 


W (I)  >  0 


1*2, 


(1.16) 
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so  that  the  Baker-Ericksen  inequality  is  not  violated.’ 

The  linear  theory  of  elastostatic  plane  strain  is  recovered  from 
the  finite  deformation  theory  briefly  described  above  by  a  systematic 
linearization  with  respect  to  the  displacement  gradients  u  Under 

Ot»P 

this  linearization,  the  distinction  between  true  and  nominal  stresses 
disappears,  and  the  constitutive  law  passes  over  into 

Ta$  =  aaB  =  2v  YaB  "  p6ctS  * 

where 

?<Ve  +  ue,a> 

are  the  components  of  the  infinitesimal  strain  tensor,  and 

U  =  2W'(2)  (1.19) 

is  the  infinitesimal  shear  modulus.  The  incompressibility  condition 
J  *  1  linearizes  to 

Yaa  *  ua,a  =  div  !*  =  0  *  (1-20) 

The  approximate  form  of  W  for  infinitesimal  deformations  is  found  by 
linearization  to  be 

W  =  £  Ya0  Ya6  ’ 

A  deformation  (1.1)  of  the  form 
’see  [7]. 


(1.17) 

(1.18) 


(1.21) 
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L 


yo  =  BctB  x$ 


0.22) 


where  the  Ba8  are  constants  satisfying 

detCBa$]  =  1  ’  0-23) 

is  a  homogeneous  deformation  of  the  incompressible  body.  Two  particu¬ 
lar  homogeneous  deformations  are  of  special  interest:  uniaxial  stress 
and  simple  shear.  For  the  former  one  takes 


yl  =  XX1 


X  >  0 


(1.24) 


with  x  constant.  From  (1.13),  one  then  finds  that  =  r ^  =  0,  and, 
if  p  is  chosen  to  be 

p  =  2W 1 ( I )  x"2  ,  (1.25) 

where 

I  =  X2  +  X”2  ,  (1.26) 

one  has 

t22  -  0  ,  (1.27) 

as  well.  The  only  nonvanishing  stress  component  is  then  found  from 
(1.13)  to  be 

xn  =  2W'(I)(X2-  X'2)  .  (1.28) 

For  simple  shear,  one  has  the  homogeneous  deformation 

t 

y*i  —  1^2  »  y2  s  *2  »  (i  .29) 
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where  the  constant  k  is  the  amount  of  shear.  From  (1.13)  one  finds  the 
relation  between  the  true  shear  stress  and  the  amount  of  shear  k 
to  be 

t12-2W(I)k,  (1.30) 

where  now 

I  =  2  +  k2  (1.31) 

We  shall  assume  throughout  that  W  has  the  following  property: 

W(I)  =  AIn  +  o(In)  as  I  «  .  (1.32) 

where  A  and  n  are  material  constants  satisfying 

A  >  0  ,  n  >  1/2  .  (1.33) 

For  an  incompressible  material  satisfying  (1.32),  one  sees  from 
(1 .26) ,(1 .28)  that  in  extreme  uniaxial  stress  (X  -*-“),  one  has 

T'j i  'v*  2n  A  X2n  ,  X  -*■  ®  .  (1.34) 

For  severe  simple  shear  (k-*-®),  one  obtains  from  (1 . 31 ) , ( 1 .30),(1 .32) 

x12  o,  2nAk2n_1  ,  k  -*■  ®  .  (1.35) 

Since  n  >  1/2,  the  stress  response  in  uniaxial  stress  is,  by 
(1.34),  always  asymptotically  hardening  as  X  -*•  ®,  in  the  sense  that 
dx^/dX  is  increasing  with  increasing  X.  In  shear,  the  stress  re¬ 
sponse  of  (T.35)  is  hardening  as  k  +  «  for  those  materials  with  n  >  1, 
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softening  for  ^  <  n  <  1,  and  asymptotically  linear  If  n=l.  The 
asymptotic  forms  of  the  stress  response  curves  in  uniaxial  stress  and 
simple  shear  for  materials  satisfying  (1.32), (1.33)  are  shown  in  Figure  1. 


If  one  were  to  permit  n  <  1/2  in  (1.32)  one  would  find  that  the 
field  equations  of  the  equilibrium  theory  would  cease  to  be  an  elliptic 
system  at  sufficiently  severe  deformations;  see  [6]. 

Before  proceeding  to  the  specific  problems  to  be  discussed,  it  is 
useful  to  take  note  of  an  implication  of  the  field  equations  (1.3), 
(1.9),  (1.14).  One  can  show  that  det  F  =  1  implies  that 


F  a  -  0  on  II 
By  ap  py  ,3 


(1.36) 


Substitution  from  (1.14)  into  the  equilibrium  equations  (1.9)  then 
gives,  with  the  help  of  (1.36),  the  equation 


[2W (I)  F  ]  R  =  P  o  e  F  on  n 
ap  ,B  ,B  By  OP  PY 


(1.37) 


If  one  multiplies  (1.37)  by  F^,  makes  use  of  the  fact  that  det  F  =  1 
as  well  as  of  the  definitions  (1 .4)^ ,(1 .5)^ ,  one  finds  that 


Vp  =  2W ' ( I )  FTV2y  +  2W"(I)  FTF  VI  on  n  . 


(1.38) 


We  will  find  this  form  of  the  equilibrium  equations  helpful  in 


the  sequel. 
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2.  THE  HALF-PLANE  DEFORMED  BY  A  CONCENTRATED  FORCE 
A.  Formulation  of  the  Problem 

We  consider  the  case  in  which  the  open  cross-section  II  of  the  un¬ 
deformed  body  is  the  half-plane  >  0,  -«  <  x2  <  ~,  and  we  denote  by 

n  the  closure  of  II  with  the  origin  deleted.  Given  the  plane  strain 
elastic  potential  W(I)  of  the  homogeneous,  isotropic,  incompressible 
material  to  be  considered,  we  seek  a  deformation  yQ  =  ya(x.|  ,X£)  on 

O 

H  such  that  the  nominal  stresses  aa6  generated  by  the  deformation 
through  (1 . 2 ) - ( 1 .5)and  (1.14)  conform  to  the  equation  of  equilibrium 
(1.9).  We  further  assume  that  the  free-surface  conditions 

<J11(0,x2)  *  a21(0,x2)  =  0  ,  |x2|  >  0  (2.1) 

hold  and  we  require  that  ,  as  Ixl  the  true  stress  field 

1  *+* 1 

should  tend  to  zero: 

Tae<xl’x2)^°  as  M  “  *  X1  >  0  .  (2.2) 

Further,  we  Impose  the  requirement  that 

oag=0(r^)  as  r+0,  uniformly  in  9,  -  <.  0  <.  ^  ,  (2.3) 

where  r,9  are  polar  coordinates  at  the  origin:  x-j  =  r  cos  0,  x2  *  r  sin  9. 
We  next  prescribe  that 
it/ 2 

{  actf}  Vd8  *  F  Sla  ’  r  >  0  (2-4) 

-it/ 2 

corresponding  to  a  concentrated  force  of  magnitude  |F|  acting  on  the 
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boundary  of  the  deformed  body  in  a  direction  parallel  to  the  x-j-axis. 

We  shall  limit  our  attention  to  the  case  F  >  0,  so  that  the  force  is 
in  the  negative  x, -direction  and  is  therefore  tensile.  In  (2.4),  n  is 
the  unit  vector  in  the  radial  direction.  The  compressive  case  could  be 
treated  by  a  similar  analysis. 

We  finally  assume  the  elastostatic  field  to  be  symmetric  about 
the  x^-axis.  This  in  particular  rules  out  a  concentrated  moment  at 
the  origin;  symmetry  also  implies  that  (2.4)  holds  automatically 
for  o=2. 


B.  The  Elastostatic  Field  near  the  Origin— Lowest  Order  Asymptotic 
Analysis 

We  now  assume  that  the  elastic  potential  W(I)  satisfies  (1.32), 
(1.33),  and  we  investigate  the  local  structure  of  the  field  near  the 
point  of  application  of  the  force.  We  begin  by  making  the  Ansatz 

m  m  , 

ya  *  r  a  va(0)  +  o(r  “"  )  as  r-*>0,  (no  sum  on  a)  ,  (2.5) 


uniformly  for  — ir/2  <.0  <_rr/2,  where  m^  and  mg  are  constants  restricted 
by 


m.|  <  1  ,  mg  >  1 


(2.6)1 


and  neither  of  the  unknown  functions  vje)  e  C2([-  J ,  *])  vanish  identi¬ 
cally.^  Moreover,  in  view  of  the  prevailing  symmetry  one  has 

One  can  show  systematically  that  (2.6)  are  necessary  when  the  applied  force 
Is  tensile.  They  correspond  physically  to  the  fact  that,  along  the  line  of 
symmetry  9*0,  the  principal  stretch  in  the  x^ -direction  is  large,  while 
Xg  is  small .  The  hypothesis  (2.6)  must  be  altered  when  the  load  is  compres¬ 
sive.  Note  that  we  do  not  assume  m^  >0. 

2 

We  actually  need  the  slightly  stronger  assumption  that  v-j ( e )  t  0  and 
Vg  has  at  most  a  finite  number  of  zeros  in  [-  £,  £]. 
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v^  (-9)  =  v^e)  ,  v2(-0)  =  -v2(0)  .  (2.7) 

It  Is  assumed  that  (2.5)  may  be  formally  differentiated  twice. 

From  (2.5)  one  obtains  for  the  deformation  gradient  tensor  F  the 

local  asymptotic  representation 

ma-l 

F  o  f  o  r  as  r  -*•  0  (no  sum  on  a)  »  (2.8) 

ap  exp 

provided 

fcx8  =  mava(9)  V6)  +  eYecY(e)  ''c^9*  ^no  Sm  on  *  (2*9) 

Here  the  dot  denotes  differentiation  with  respect  to  0  and  we  have 
introduced  the  abbreviations 


c^ (9)  =  cos  0  ,  c2(0)  3  sin  0  . 


(2.10) 


From  (2. 8), (2. 9)  there  follows 


J  s  det  F  =  (mjV^  -  m2v2Vj)  r 


m.|  +  m2  -  2 


m,  -Htu  -  2 

+  o(r  c  )  ,  r  -*•  0 


(2.11) 


Since  incompressibility  requires  J  =  1,  we  must  have 


m-j  +  m2  -  2  <_  0  . 


(2.12) 


and  either 


m^v^v2  -  m2v2v-|  =0  if  m1+m9  <  2  , 


1  "’2 


(2.13) 


or 


MfT|  ■,[  .  BUM  mmUmm 
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ml v1^2  "  m2v2^1  *  1  if  ml+m2  =  2 


(2.14) 


From  (1.5), (2. 8), (2. 9)  we  obtain 

2(mrl) 


hr 


where 


6(0)  as  r  -*■  0 


G(0 )  *  v^(8)  +  v^(e) 

In  view  of  the  assumption  (2.6)  concerning  m^,  one  has  from  (2.15)  that 
I  -+■  od  as  r  -*■  0.  The  material  assumption  (1.32)  then  yields 


(2.15) 


(2.16) 


2n(m.-l) 

W(I)  *  AGn(0)  r 


W'(I)  *  nAGn  (6)  r 


.  2(n-l)(m1-l) 

n-l/«\  _  I 


y  as  r  0 


W"(I)  *  n(n-l)  AGn  (0)  r 


2(n-2)(m.-l ) 


n-2,^  _fc'"  w/v"l 


V 


(2.17) 


We  now  recall  the  field  equations  in  the  form  (1.38);  with  the  help  of 
(2. 17), (2. 8), (2. 9), (2. 6), (2. 15)  we  find  from  (1.38)  that 

2(m-j-l  )n  -  1 


2nAm^v^(9)  Z(9)  r 


and 


iff  *  2nAv1(0)  Z(9)  r 


2(mrl)n-l 


(2.18) 


(2.19) 


as  r  -►  0  ,  where 
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Z(8)  =  6n-2(0)  {G(0)[v1(0)  +  mf  Vl(0)] 

+  (n-l)[G(0)  v,(0)  +  2m.| (n-1 )  G(0)  v^©)]}  .  (2.20) 

Compatibility  of  (2.18), (2. 19)  requires  that  and  Z  satisfy 

m-jVjlt  +  [1  -  (2n-l ) (m^-1 )]  v-jZ  =  0  .  (2.21) 

Once  (2.21)  is  fulfilled,  one  finds  from  either  (2.18)  or  (2.19)  that 
m,  2n(m,-l) 

p  •'*  A  v.(8)  Z(0)  r  as  r  -*■  0  •  (2.22) 

lTl«j  “  I  I 

We  next  consider  the  boundary  conditions  (2.1).  Because  of  (1.14) 
these  are 

2W(I)F11  -  pF22  =  0 

•  at  0  =  ,  r  >  0  .  (2.23) 

2W’(I)F21  +  pF]2  =  0  J 

Multiplying  the  first  of  (2.23)  by  F^,  the  second  by  F^,  adding  the 
results,  and  using  (1.3),  we  obtain 

p  =  2W,(I)(F21  +  F21)  at  0  -  ±|,  r  >  0  •  (2.24) 

On  the  other  hand,  eliminating  p  between  the  two  equations  (2.23)  yields, 
in  view  of  (1.16), 

^11^12  ^  ^21^22  *  ®  3t  0  *  i ^ i  r  >  0  •  (2.25) 


Making  use  of  (2. 17), (2. 8),  and  (2.9),  we  obtain  from  (2.24)  the  result 
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„  i  2n(m.-l ) 

p  2n A  Gn"  (0)  (0)  r  as  r  -►  0,  0  *  ±  j  *  (2.26) 

Comparing  (2.22)  at  0  *  ±j  and  (2.26)  leads  to 

2nGn“1(±f-);i(±f)  3  ijTT  Vl(±|)Z(±J)  .  (2.27) 

The  second  boundary  condition  (2.25),  with  the  help  of  (2. 8), (2. 9), (2. 15) 
gives 

¥,(*})  )  -  0  .  (2.28) 

We  now  show  that  the  two  boundary  conditions  (2.27)  and  (2.28), 
together  with  the  differential  equation  (2.21),  imply  that 

v^if)  *  0  (2.29) 

and 

Z(±  |)  =  0  .  (2.30) 

We  first  establish  (2.30)  by  showing  that  the  hypothesis  Z(ir/2)  f  0  leads 
to  a  contradiction.  If  Z(tt/2)  f  0,  there  is  an  interval  [0o,tt/2]  , 

0O  <  ff/2,  on  which  Z(0)  vanishes  nowhere,  by  continuity.  Equation  (2.21) 
can  then  be  integrated  on  [0o,tt/2]  to  give 

ml 

Vj(9)  *  C  Z(0)  1  ,  0oi9<J»  (2.31) 

where  C  Is  a  constant.  If  C  *  0,  then  v^  (9 )  =  0  on  [0q,it/2],  whence  by 
(2. 16), (2.20),  Z(0)  =  0  on  [90,ir/2],  contradicting  the  hypothesis 
Z(tr/2)  f  0.  Thus  Z  f  0,  and  so  by  (2.31),  v^O)  f  0  for  all 
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0  e  [9q.1t/2].  In  particular,  v^(ir/2)  t  0,  so  that  by  (2.28),  (tt/2 )  =  0 . 

Thus  (2.27)  has  been  violated  unless  m.|  =  0.  But  if  =  0  we  conclude 

from  (2.21)  that  =  constant  on  [0q,it/2],  which  leads  via  (2.20) ,(2.16) 
to  Z  =  0  on  [90,ir/2],  again  contradicting  the  hypothesis.  Thus  indeed, 
Z(ir/2)  =  0,  and,  since  Z(8)  is  even,  (2.30)  holds.  From  (2.27)  and 
(2.16)  It  then  follows  that  (2.29)  holds  as  well. 

An  argument  similar  to  that  just  used  to  establish  (2.30)  can  now 
be  constructed  to  show  that  (2. 21), (2. 30)  imply  that 

Z(0)  =  0  on  [-it/2,tt/2]  (2.32) 

From  (2. 32), (2. 20), (2. 29)  we  then  obtain  a  nonlinear  eigenvalue  problem 
for  m^ .Vj (0): 

[Gn*^ (9)  ^(8)]*+  [m^  +  2m1(m1-l)(n-l)]Gn"1(e)  v^e)  =  0  , 

-ir/2<e<ir/2  (2.33) 

v1(±  f)  =  0  .  (2.34) 

The  differential  equation  (2.33)  is  identical  with  one  which  has  arisen  in 
the  local  analysis  of  the  elastostostatic  field  near  the  tip  of  a  crack; 
see  [8], [9], [10],  [11] 

In  view  of  (2.32),  we  have  from  (2.22)  that 

2n(m,-l ) 

p  *  o(r  )  as  r  -*•  0  .  (2.35) 

Suppose  that  m^  *  0.  Then  (2. 33), (2. 34)  imply  that  v^(9)=  constant 

m-i 

and  the  leading  term  r  v^(8)  in  the  expansion  of  y^  near  r  *  0  may  be 
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vlewed  as  a  rigid  body  translation  parallel  to  the  x^-axis.  Since  the 
boundary  value  problem  determines  the  elastostatic  field  at  best  to 
within  an  arbitrary  translation  of  this  kind,  we  shall  discard  the  case 
m^  =  0  and  assume  henceforth  that,  in  addition  to  (2.6), 

m-|  f  0  ,  (2.36) 

also  holds. 

It  is  now  possible  to  prove  that  (2.36), (2.6)  and  the  assumptions 
made  concerning  v^  and  imply  that  (2.13)  leads  to  the  contradiction 
v-|  =  0.  Thus  (2.14)  must  hold,  and  thus 

m^  +  -  2  .  (2.37) 

From  (1.14)  we  have 

0i8  =  2W(I)F16-pe6YF2Y.  (2.38) 

Making  use  of  (2.17),(2.8),(2.9),(2.35)-(2.38)  we  can  show  that  the  first 
term  on  the  right  in  (2.38)  dominates  the  second,  and  hence  that 

n-i  (2n-l )(m,-l ) 

Ojg  ^  2nAm1Gn  (9)  f^(8)  r  as  r  -*•  0  .  (2.39) 

From  (2.3)  we  conclude  that  (2n-l)(m^-l)  >.  -1.  In  order  to  use  (2.4) 
with  o  *  1,  we  first  observe  that  ng  *  c^  (see  (2.10)),  and  from  (2.39), 
(2.9)  that 


I 
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I/2  yZ  n-l  (2n-l)(m.-l)*l 

j  aigng  r  d0  ^  2nm^A  G  (0)  v^(0)  d0  r 

'v/Z  "ir/2  as  r  -*•  0  .  (2.40) 


It  follows  that 

rol  s  ^TT^*  <  1  n  M  .  (2.41 11 

From  the  boundary  value  problem  (2. 33), (2. 34)  with  m-j  given  by  (2.41),  one 
finds 

v-j(0)  *  *  constant  ,  -ir/2  <.0  <_  tt/2  ,  n  f  1  .  (2.42) 

Since  (2.37)  holds,  we  have 

m2  =  2’ral  =  2n?T  >  1  ’  n  f  1  .  (2.43) 

With  the  help  of  (2.42)  and  the  fact  that  Vg(0)  is  odd,  we  can  now  deter¬ 
mine  Vg  from  (2.14)  as 

v,(e)  *  TT“  8  =  1^9  9  »  -ir/2  £  0  £  Tt/2  ,  n  f  1  .  (2.44) 

i  K]m1  * 

Finally,  we  return  to  (2.4)  with  a  =  1  to  determine  in  terms  of  F. 

Using  (2 . 40 )- (2 . 42 )  and  (2.16),  we  obtain 

r(miki)2]  (miki>  *  » 


^Recall  that  m1  a  0  has  been  excluded. 


n  t  1 


(2.45) 


p 


and  hence 
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k,  = 


2n-l 


1  '  2(n-l )  '2n  Air 


f  O  n  flw) 


2n-l 


n  f  1 


(2.46) 


Thus,  for  n  t  1,  we  have  determined  the 
proximation  to  the  deformation  near  r  =  0  as 


y,  •v  k,  >• 


2(n-!l 

2n-l 


first  terms  (2.5)  in  the  ap- 
follows: 


y2 


2n-l  1  __ 

k7r 


2n 


y  as  r-*0,  -ir/2  <.  8  <.  tt/2  ,  n  f  1  , 

(2.47) 


with  k1  related  to  F  through  (2.46).  The  deformation  image  of  the  bound 
ary  0  =  ±ir/2  of  the  half-plane  is  then  given  in  first  approximation  by 


yl 


1  fJ_iT/n  i„  |1  -i 

ir  (4nA}  ly2'  n 


|y2l  -  0  .  n  M  .  (2.48) 


We  note  that  if  n  >  1,  so  that  the  material  is  asymptotically  hardening 
in  simple  shear  [see  (1.35)  and  Figure  1],  the  displacement  under  the 
load  is  finite,  while  this  is  not  the  case  for  softening  materials  (n  <  1 
A  sketch  of  the  deformed  surface  based  on  (2.48)  is  given  in  Figure  2. 

The  case  n  *  1  (a  material  which  is  asymptotically  linear  in  shear 
(Figure  1))  has  been  excluded  in  the  results  (2. 47), (2. 48).  To  treat 
this  case,  it  is  necessary  to  replace  the  Ansatz  (2.5),  (2.6)  by 

y-|  *  (log  r)  v^Q) 
nu 

y2  *  r  v2(0)  ,  m2  >  1 


as  r  -*■  0  . 


(2.49) 
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The  special  nature  of  the  case  n  =  l  arises  because  =  0 
and  =  2(n-l  )/(2n-l  )"are  both  eigenvalues  (adjacent  ones  in  fact)  of 
the  problem  (2. 33), (2. 34).  They  are  distinct  as  long  as  n  f  1 ,  but 
coalesce1  as  n  +  1 .  This  coalescence  may  be  used  to  motivate  the  form 
of  the  new  Ansatz  (2.49)  for  n  =  1;  we  omit  the  details.  One  finds  from 
(2.49)  that  m^  *  2,  v-j(e)  =  k1  =  constant,  Vgte)  =  (1/^)8  ,  where 
k.|  =  F/2Air.  The  counterparts  of  (2.47)  are 

r  ' 

*i  %  m log  r 

r  as  r  -*■  0,  -ir/2  <.  0  <.  tt/2  ,  n  =  1 , 
y2  ^  9  »  I  (2.50) 

while  the  deformed  boundary  is  now  given  approximately  by 


*1 


log( 


as  |yg  1  -*>  0  ,  n  =  1  . 


(2.51) 


The  displacement  is  unbounded  near  the  point  of  application  of  the  load, 
as  It  is  for  the  softening  material  (n  <  1). 

Although  the  nominal  stresses  0^,0^  are  fully  determined  to  lead¬ 
ing  order  as  r  0  at  this  stage,  the  fact  that,  as  yet,  only  the  weak 
estimate  (2.35)  is  available  for  the  hydrostatic  pressure  p  makes  the 
asymptotic  determination  of  CT22 *CT21  ’mP0SS1^e  without  higher  order 
considerations.  In  view  of  the  relationship  (1.7)  between  the  nominal 
stresses  0Qg  and  the  true  stresses  Tag,  the  full  asymptotic  determination 

1 A  similar  but  more  complicated  coalescence  of  eigenvalues  arises  in 
crack  problems,  see  [8],  [9]. 
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of  the  latter  must  also  await  such  higher  order  results. 


C.  Higher-Order  Asymptotic  Considerations 

For  the  present,  suppose  that  n  /  1  and  replace  (2.5)  by  the  two- 
term  asymptotic  representations 
m,  s. 


1  51  *1 
y-j  ^  k +  w^ (6)r 


s„ 

y2  %  k20r  +  w2(9)r  c  ,  J 


>  as  r  -*•  0  ,  -it/ 2  <  0  i  tt/2  ,  (2.52) 


with  the  stipulation  that 

5 1  ^  )  02  ' 


s-i  >  m,  ,  s,  >  m9  ,  wae  C2([-|-,  |])  ,  wa  jf  0  on  |]  , 


(2.53) 

and  that  w^ ,w2  have  the  respective  parity  of  v1  and  v2>  Equations  (1.2) 
(2.52)  lead  to  the  following  representation  for  the  components  of  the 
deformation  gradient  tensor: 

m i  -l  s  -1  s  -1 

Fae  s  fa6  r  +9a0  r  *  +°(r  0  )  »  (™>  sum  on  a),  (2.54) 


provided  fag  is  given  by  (2.9)  and 


9o6(0)  =  sacg(9)  wa(e)  +  e  gcy(e)  wa(0)  (no  sum  on  a).  (2.55) 

The  asymptotic  representation  of  the  deformation  invariant  I  depends  on 
the  value  of  s-j .  One  can  show  after  some  calculation  that  necessarily 

m.|  <  s1  <  4-  3m.j  .  (2.56) 
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Then  the  asymptotic  representation  for  I  is 

2  2  )  m.  +  s,  -  2 

1 '  "i  »i r  +  2misiviv,i r  • 

Equations  (2. 17), (2. 54), (2. 55), (2. 9), (2. 57)  and  (1.38)  yield 


s,-3 


§£*  2nAGf_1  Y(e)  r  1 

l2£-0(rSl'3) 


r  39 


r  r  -*•  0  , 


where 


and 


Y(e)  =  +  kw^  , 

<  *  s1C(2n-l)s1  -  2(n-l)]  , 

Q-|  (0 )  *  MjVjte)  =  m1k1  . 

On  the  other  hand,  the  boundary  conditions  (2.23)  lead  to 

2(s,-l )  -  m, 

W1  (*  f)  =  0  *  P(r,  ±|)  ^  o(r  1  ])  . 


(2.57) 


(2.58) 


(2.59) 


(2.60) 


(2.61) 


Integrate  the  first  of  (2.58)  with  respect  to  r  to  get 
2n  A  G?n_1  s,-2 

- Y(9)  r  ,  (Sj  f  2.)  .  (2.62) 

* 

Comparing  (2.62)  with  the  second  of  (2.61),  one  deduces  that 

Y(±|)  *  0  .  (2.63) 


The  compatability  of  (2.58) ,(  2.62)  together  with  (2.63)  then  gives 
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Y(0)  =  0  on  [-ir/2,ir/2]  .  (2.64) 

Equations  (2.59)  and  the  first  of  (2.61)  imply 

w-j  (9 )  =  B  cos  Sk  9  ,  B  constant,  -  ir/2  <9  <  tt/2  ,  (2.65) 

Sk  f  =  jif  »  j  =  0,  ±1 ,  ±2,*  •  •  .  (2.66) 

The  second  of  (2.59)  and  (2.66)  imply  that 


,  -  2(n-1)  t  J[2(n-l)f  +  16(2n-l)  j2  .  . 

1  2(2n-l ) 

One  seeks  the  smallest  value  of  s1  satisfying  (2.56). 
j  =  1 ,  so  that 


±1,  ±2,  •••  .  (2.67) 

This  occurs  for 


S1  =  2TTT +  7^7^  *  (2’68) 

The  asymptotic  results  for  the  spatial  coordinates  deduced  this  far 
may  be  summarized  as  follows: 

gfcU  s 

y^  ^  r  2n”  +  B  cos  20  r  1  ,  n  /  1  ,  (2.69) 

2n 

y2^k2e/n"'  ,  nM  ,  (2.70) 

sr2 

p  ^  o(r  '  )  ,  n  >  1/2  .  (2.71) 

The  case  n  =  1  is  treated  at  the  end.  Substituting  from  (2.54)  into  (1.3), 
using  (2. 11), (2. 14)  gives  that 


s,  +m2  -  2 


1 

J  'v  1  +  (k2siwi  “  mgkgW^Q)  r 

s,  +  m,  -  2  s,  +  s,  -  t 
+  m.jk.|W2  r  2  1  +  0(r  1  2  )  . 


(2.72) 


I 

'  1 1  ‘MTittWUteBimill' 
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A  simple  analysis  gives 


s2  =  sT  +  m2  -  m1  =  s1  +  2n-1 


and 


*  =  ?2  4  0  •  .  Vl 

w2  m]  Tc^  9  W1  Fjm^  W1 


which  on  integration,  using  (2.65),  yields 


(2.73) 


(2.74) 


O  I 

w2  = - j  (m2e  cos  2Q-j  (m2  +  s^)  sin  26}  .  (2.75) 

(m-jk^ ) 

For  n  =  1,  we  can  similarly  show  that 

2 

y -|  'v  k1  log  r  +  B  cos  29  r  ,  n  =  1  .  (2.76) 

In  an  effort  to  find  a  strong  estimate  for  the  pressure  field,  we  now  as¬ 
sume  the  following  three-term  asymptotic  representation  for  the  deforma¬ 
tion: 

m  s  t 

ya  *  var  °  *  wa  r  a  +  za  r  a  ,  (no  sum  on  a),  (2.77) 

with  the  stipulation  that 

t|  >  s1  >  m1  ,  t2  >  s2  >  m2  ,  (2.78) 

and  z^,z2  are  functions  possessing  derivatives  of  second  order  on 
[-ir/2,ir/2],  which  fail  to  vanish  identically  and  have  the  same  parity  as 
v.j  and  v2. 

For  n  >  1,  It  can  be  shown  that 

Zj(8)  *  D  cos  2/t^  9  +  k|(pj92  +  u2)  ,  (2.79) 
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and 


t1  =  2m2  -  m1  = 


where 


(m2 +  2y, )  2  _ 

D  =  —*= — 1  -1  Z  ,  1  <  /t7  <  2  , 

2/ty  sin  /ty  ir  1 


u-j  =  E(2  —  3m.| )  (nig  *  )  *  4m2]  m2  * 


u2  3  -  2(m.|  +  )] 


V  *  (2-3m1)[l  -  (m2-m,)  (f)2]  ,  J 


and  we  arrive  at  a  strong  estimate  for  the  pressure  field, 
p-  2nAG^<"-1)  ^  [92 -  (f)2]  +  1 )  r  ^  ,  n 


For  n  <  1 ,  we  can  show  that 


z-, (e )  =  B2(u1  cos  40  -  y2) 


and 


where 


t1  *  2s1  -  m1  , 


y-j  3 “ ^  (2n-l )t^s^k2  ,  y2  *  ^  (n-l)s^k2  , 


but  for  the  pressure  field  we  have  only  the  weak  estimate 

2s,  -  m,  -  2 
p  *  o(r  ) 


(2.80) 

(2.81) 

1  .  (2.82) 

(2.83) 

(2.84) 

(2.85) 

(2.86) 
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At  this  stage  we  can  determine  z2,t2  ^rom  condition  (1.3);  however, 
it  Is  not  necessary  to  record  the  results  here. 

To  find  a  strong  estimate  for  the  pressure  field  for  n  <  1  we  assume 
the  deformation  admits  the  representation 
m is  t  i 

ya * v  *  v a  +  v a  +  V °  *  (no  sum on  a)’  (2*87) 

with  the  stipulation  that 

>  tj  >  Sj  >  idj  ,  SL^  >  *2  >  s2  >  m2  * 

qa  g  C2([-n/2,Tr/2])  ,  qQ  =  0  on  [-ir/2,ir/2]  . 

We  can  now  show  that 

Aj  =  2m2  -  m1  ,  7/12  <  n  <  1  ,  (2.88) 

and  that  q-j(e)  Is  given  by  the  value  of  z^e)  in  (2.79).  Thus,  we  can  see 
a  trade  in  dominance  between  the  third  and  fourth  term  of  (2.87),  for  a  *  1, 
as  n  passes  through  n  =  1.  Condition  (1.3)  is  again  used  to  determine  *2 
and  q2.  For  7/12  <  n  <  1/2,  equation  (2.82)  is  found  to  give  a  strong 
estimate  for  the  pressure  field.  The  value  n  *  7/12  is  a  transition  point 
for  the  pressure  field.  A  strong  estimate  for  n  in  the  range  (1/2,7/12) 
requires  much  further  analysis. 

At  this  point,  we  will  record  the  results  for  n  *  1. 

y^  'v  k1logr+  Bcos  29  r2  +  J-  k2cos  49  r4log  r 

k2 

+  [Ecos  40  -  (36  sin40  +  462  +  ir2-|)]r4  , 


1 
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y2  ^  k2e  r 

p  -v  Ak2(402  -  it2  +  2)  r2 


where  B  and  E  are  constants. 


(2.89) 


D.  Summary  of  Results  for  Deformation  and  Stresses 


The  asymptotic  results  for  the  spatial  coordinates  and  pressure  field 
may  be  summarized  as  follows: 


2( n-1 )  _  o„  _ 

— * o  ^  A  lS ^  -fH-.  i  y 

yi  k1  r  ^n"  +  B  cos  29  r  +  B  (y-,Cos  40  -  y2)r  ,  n<.^-  , 


^  s,  7  A  2s, -m, 

y-|  ^  k1  r  ^n“  +  B  cos  20  r  +  B^(y-,cos  40  -  y2)r 


*1  '  kl 


-  o  4-3m,  y 

+  [D  cos  2/1^9  +  k^y^+y^r  ,  n <  1  » 

s,  o  p  4- 3m, 

+  B  cos  20  r  +  [Dcos  Z/t^Q  +  k|(y.|0  +y2)]r  ,  n  >  1  , 


2n 

y2  *  k20  r2""1  +  w2(0)  r 


S1  +  27T 


,  n  f  1 


p  *  £  k^  { — [02-  (J)2]  +  1}  r  2n_1  ,  X<n<»,  n  f  1  , 


(2n-l) 


(2.90) 


where  k^.k,,  are  given  in  (2.44) ,(2.46) ;  w2(0)  in  (2.75);  s^  and  t.j  in 
(2. 68), (2. 80);  D,  y^  u2  in  (2.81);  and  yj.y2  ’n  (2.85).  B,  a  constant,  is 
left  undetermined  by  the  local  analysis. 


Me  now  turn  to  the  asymptotic  determination  of  the  actual  stresses 
From  (1.13), (2. 17), (2. 8), (2. 9), (2. 41)- (2. 44),  and  (2.35),  one  finds 


Tn'7<miki>^  • 


-  2 ( n- 1 


t12  =  T21  *  ir  ^2*2*  r 


n  f  1 


e,  n  f  1 


(2.91} 


(2.92) 


and  using  (2.82  ), 


t22  *  ff  m2k2  +  r 


2rvT  '2' 


,  ^-<  n<®  ,  n  f  1  .  (2.93) 


The  true  stresses  here  are  referred  to  the  material  polar  coordinates 
(r,0).  The  stress  component  t,.j,  which  is  of  primary  physical  interest, 
becomes  unbounded  at  the  origin  for  all  admissible  values  of  the  hardening 
parameter,  the  singularity  becoming  more  severe  with  decreasing  values  of 
n,  and  for  the  range  of  n  under  consideration  is  always  stronger  than  that 
predicted  by  the  linear  theory.  The  other  normal  stress  component  t22  're¬ 
mains  bounded  for  n  <.  2,  but  for  n  >  2  it  becomes  unbounded,  while  the 
actual  shearing  stress  t-|2  is  bounded  for  n  <  1,  and  becomes  unbounded  at 
the  origin  for  hardening  materials  (n  >  1).  In  both  stresses  t22,t-|2  .  the 
severity  of  the  singularity  increases  with  increasing  values  of  n;  however, 
the  order  of  the  singularity  is  less  than  that  predicted  by  the  linear 
theory  for  all  allowable  n. 

For  n  *  1  we  have.  In  sunmary,  that 

2  3  3  4 

y1  ^  k^log  r  +  Bcos  29  r  +  ^  kjj cos  40  r  log r 

k3 

rc  An  J.  2  tea  .  Jm  A  a  An^-  j.  _2  5  \  T 


+  [Ecos  49  +  tI-  (60  sin  40  +  49  +tt  -  5-)]  , 


y2  'v  k-0r2  +  2k2  B  (0  cos  0  -  si n  0 )  r* 


(2.94) 
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p  -V  Ak^e2  -  *2  +  2)  r2 

Here,  k^  *  l/kg  =  ,  while  B  and  E  are  constant,  undetermined  by  the 

local  analysis.  The  components  of  the  actual  stress  tensor  are,  using 
(1.13),  (2.17),  (2.8),  (2.9),  (2.50),  and  (2.35), 

Tll^klr‘2  *  -22  -^2k|{e2  +  (5-)2}  r2  , 

(2.95) 

t12  *  t21  *  7  *2k2*  9  ‘ 


At  this  stage,  a  comparison  of  the  results  given  in  (2. 94), (2. 95) 
with  those  predicted  by  the  classical  linear  theory  reveals  the  presence 
of  a  ljgr  singularity  in  the  dominant  term  of  the  deformation  in  the 
direction  in  both  cases.  This  is  the  only  similarity!  According  to 
linear  theory  all  the  components  of  the  stress  tensor  possess  a  1/r  singu¬ 
larity  at  the  origin,  but  from  (2.95)  we  see  the  stress  component 
is  more  singular,  while  the  other  components  are,  in  fact,  bounded  there. 
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